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An equilibrium theory of classical fluids based on the space distribution among the particles is 
derived in the framework of the energy minimization method. This study is motivated by current 
difficulties of evaluation of optical properties in atmospheres of degenerate stars. Present paper 
focuses on diluted one-component systems, where the interaction energy is calculated as a sum of 
binary contributions. The spatial configuration of the gas is described in terms of a particle-state 
variable v which roughly measures the free space surrounding each particle. The formalism offers 
a unified treatment of both thermodynamics and structure of fiuids, since it not only provides the 
state function (the Helmholtz free energy) of a fluid, but also automatically gives information on 
the microstructure of the system (e.g. the nearest-neighbor distribution function). Equations of 
state and nearest-neighbor distribution functions of perfect and hard-body fiuids are obtained in 
straightforward way in one, two and three dimensions. The formalism is applied to describe the 
atomic population in a partially ionized hydrogen gas. Combined and self-consistent evaluations 
of atomic populations and internal effects on bound states are performed in a detailed form. The 
present theory allow us to resolve the atomic density at each internal level into groups of atoms 
which experiment different perturbation intensities according to the size of their spaces v. 
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I. INTRODUCTION 

The evaluation of the thcrmodynamical, transport and 
optical properties of non-ideal plasmas is of great inter- 
est to characterize the behavior of the matter at different 
physical conditions. It has many applications in mate- 
rial sciences, geophysics and astrophysics. In particular, 
the equations of state provide the basic thermodynamic 
quantities needed to determine the physical properties 
of stellar envelopes and interiors. Models of gases must 
also provide detailed atomic and molecular populations 
required to obtain the monochromatic opacity, a quan- 
tity which is essential to calculate realistic spectra emit- 
ted from astrophysical objects. Our study is motivated 
by current difficulties of evaluation of gas opacities in 
envelopes of white dwarf stars (WDs). 

While non- ideal effects on equations of state (e.g., 
departures from the ideal P-T-p relationship) become 
present at densities near to ionization pressure (which 
occurs above p ~ 0.3 g cm~'^ for a pure hydrogen plasma 
0), signatures of non-ideal effects (i.e., particle interac- 
tions) are evident in spectra emitted by gases at very 
low densities. As it is well-known since the earliest days 
of the quantum theory, interparticle perturbations are 
responsable for the reduction of the ionization energies, 
shifts of lines and modifications of their profiles "2^ . 

The reduction of the ionization potential along with 
the line merging, yields the advance of the series limits 
in the spectrum and forms 'pseudocontinuums', which re- 
semble extensions of the bound-free adsorptions toward 
longer wavelength Sj- Pseudocontinuum opacities may 
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produce strong modifications in the spectra of cool WDs 
with rich- hydrogen envelopes However, so far, the 

model atmospheres fail to reproduce in detail the photo- 
metric observations of WDs with effective temperatures 
below 4000 K, and the source of discrepancies likely lies 
in the computation of gas opacities which are affected by 
non-ideal effects iQ]. 

Current optical simulations 0, EL commonly used in 
the evaluation of WD spectra (e.g. P,'^,'^,'!^), are based 
on the occupation probability formalism elaborated by 
Hummer and Mihalas (HM) 11]. On the basis of ear- 
lier works and the implementation of the technique 
of free energy minimization, the HM formalism assigns 
to each atomic state j a probability Wj of finding the 
atom in this state relative to that of finding it in a sim- 
ilar ensemble of non-iteracting particles. For the com- 
putation of optical properties, the quantity Wj is consid- 
ered as an estimate of the number of states of type j 
that are available to be occupied {bound level), whereas 
1 — Wj is a measure of the fraction of j-states that are 
severely perturbed by plasma interactions, such that an 
atom in this fraction is actually unbounded (dissolved 
level) 0- In this heuristic scheme, a radiative excitation 
of an atom from a bound level to a dissolved level con- 
stitutes a pseudo-continuum absorption. Althought this 
proposal seems intuitively right, we have recently showed 
that it is not self-consistent with the HM formalism 
Briefly, the optical simulations assume the existence of 
a fraction of atoms with perturbed electronic energies 
(values shifted relative to the continuum), whereas the 
occupation HM formalism has been developed for par- 
ticles species which have unperturbed energies (ioniza- 
tion potential of isolated atom). As a consequence, the 
atomic population which contributes to a pseudocontin- 
uum absorption is overvaluated, because the abundances 
of perturbed atoms are calculated with lower bound ener- 
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gies than those assumed in their transitions to dissolved 
states. This may explain why such optical simulations 
yield a huge (unphysical) pseudocontinuum-Lyman opac- 
ity for cool WDs Q. Since the starting level of these 
processes is the fundamental state, the overevaluation of 
the Lyman pseudocontinuum is more severe than those 
of Balmer and Paschen. 

While a group of gas properties (equations of state 
and response functions, e.g., specific heats) may be cal- 
culated directly from the global distribution of particles 
(mean densities), others, such as emissivities and opaci- 
ties in general depend on particular population subsets. 
As an appropriate example, we mention the point of view 
adopted by the statistical theory of line broadening by 
pressure effects. There, the spectrum of a broadened line 
is determinated from a statistical sum, where the contri- 
bution of each radiating atom depends on the configura- 
tion of its surrounding particles . This picture can be 
also applied to continuum transitions. Thus, an appro- 
priate treatment of pseudo-continnum opacities could be 
developed on the basis of a formalism able to distinguish 
groups of atoms which, at a given instant, experiment 
different degrees of interparticle perturbation. 

Motivated by these considerations, we believe it is 
physically significant and useful to develop a statistical 
mechanical formalism for gases based on the distribution 
of the physical space among the particles. The central 
idea is that to each particle is assigned a space region 
V nearly related with the closeness of particle neighbors, 
and which is treated as a configurational parameter in 
the thermostatistical description of the gas. A particle 
with small space v will have very near neighbors and, 
therefore, it will be subjected on the average to high 
perturbations, while a particle will be roughly isolated 
if it has a huge space v assigned. This description of- 
fers an appropriate framework for the study of optical 
properties of gases. So far, the application of statistical 
mechanics in the analysis of space partition into a set of 
volumes has been circumscribed to some particular sys- 
tems such as cristals Similar theory for gases has 
not been formulated as far as we are aware. As a first 
step, the formalism developed in this paper is devoted to 
one-component gases at the limit of low densities. 

We formulate the theory in terms of the widely used 
free-energy minimization method |l5j |. Our aim is there- 
fore to develop a Helmholtz free-energy model for gases 
which contains a notion of space per particle, and to find 
the statistical equilibrium state and the thermodynamic 
properties of the system from the energy minimization 
with respect to the particle number density, subjected 
to specific conditions. Since the present work concerns 
gases at low densities, we assume factorization of the free 
energy into configurational, translational and particle- 
interaction contributions. 

The paper is organized as follows. Sec. ITU brieflv re- 
views the importance of spatial statistical in numerous re- 
search areas, and focuses on the application of space par- 
tition schemes to the analysis of spatial micro-structures 



in gases. In Sec. IIIII we introduce the volume distribu- 
tion among particles as a new configurational parameter 
in thermostatistics of gases. The gas is treated with the 
Maxwell-Boltzmann statistics and the particle interac- 
tions are represented by additive pair potentials. The 
equilibrium states are derived in Sec. IIVI In Sees. IVIVII 
the theory is applied to perfect and hard-sphere gases. 
We show that they have the expected thermodynamic 
properties. There, the Van der Waals model is introduced 
as an illustration of some of the considerations underlying 
this article. In Sec. I VIII we explore the application of the 
theory to the description of atoms perturbed by plasma 
effects in a partially ionized, hydrogen gas. Concluding 
remarks are given in Sec. IVIIII 



II. SPATIAL STATISTICS 

One of the most basic properties of a collection of ob- 
jects is the division of space among them. The partition 
of space is commonly associated to processes that create 
random geometries. There is a considerable variety of 
physical systems which show stochastic structures. For 
example: biological tissues |l6|: landscape geometries in 
agriculture and forestry !l7| ; territoriality zoology [l8l| ; 
grains in foams and metallurgial aggregates (IJ, Jl^ ; in - 
termolecular cavities in fluids [IJl; solar granulation |2ll |: 
fragmentation of interestellar clouds 22]; stellar popula- 
tions |23| ; matter and void distributions in the large-scale 
structure of the Universe [l^ . Most of these systems have 
often been studied with computer simulation of stochas- 
tic, space-filling structures. A variety of geometric con- 
structions have been applied in these analyses |0| . 

Gases in thermodynamic equilibrium are also typical 
disordered systems. In fact, althougth macroscopically 
homogeneous in the absence of external fields, these sys- 
tems have spatial inhomogeneities at microscopic levels. 
The particles of a fluid move continuously and the instan- 
taneous spatial concentration of material is non-uniform 
and fluctuating at different scales. Fig. ^ (above on the 
left) shows a group of objects random distributed in a 
surface, which represents a typical configuration of a two- 
dimensional ideal gas at a certain instant. We note that 
some particles are close enough among them while oth- 
ers are separated from their nearest neighbors by larger 
distances. One can then imagine a division of the space 
into a set of regions, one region per particle taking into 
account the closeness of its neighbors. The most likely 
well-known partition of space is given by the Voronoi tes- 
selation '2^. The Voronoi region of a given object is the 
collection of points in space nearer to that object than to 
any other. The complete set of Voronoi regions form the 
so-called Voronoi diagram of the system. An example is 
given in Fig. ^ (above on the right). There are, however, 
many other types of space divisions (e.g., generalizations 
of the Voronoi tesselations [13 )• 

It should be also observed that while the particle dis- 
tribution of a perfect gas corresponds to a Poisson dis- 
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FIG. 1: Above: Poisson distribution of 400 points in a plane 
(on the left) and the Voronoi tessellation corresponding to 
points within the dotted box (on the right). Solid lines denote 
the boundaries of the region of space closest to each particle. 
Below: Typical configuration of centers of 400 hard discs for 
a packing fraction (fraction of surface covered by the discs) 
ri = 0.2 (on the left) and rj = 0.5 (on the right). 



tribution of points in the space, gases formed by inter- 
acting particles are characterized by other different dis- 
tributions. Thus, for example, the well-known one com- 
ponent plasma (model of discrete positive charges in a 
continuous negatively charged background '2^) shows a 
long-range ordered (or 'superhomogenous') arrangement 
of particles, which can be appreciated from the analysis 
of the so-called power spectrum of the particle distribu- 
tion . Also the hard-body fluids show particular spa - 
tial patterns which have been studied extensively [3Cll| . 
As an example, the lower panels in Fig. ^show typical 
configurations of hard-disc fluids, ft is evident that the 
density fluctuations of a gas in equilibrium depend on the 
type of interactions among the particles. Furthermore, it 
is expected that the microscopical spatial structure of a 
gas in equilibrium corresponds to the most probable one 
according to the principles of the statistical mechanics. 

The present work is based on the assumption that the 
microscopical structure of a gas can be described by some 
class of partition of space among the particles, and that 
the most probable partition can be inferred using the for- 
malism of the classical statistical mechanics. The region 
assigned to a particle, denoted v, will be called the avail- 
able volume of the particle because it gives a notion of 
the size of the particle- free space which surrounds it [U ■ 
We have mentioned that the division of the space occu- 
pied by the gas is not unique and, therefore, an exact 



definition of v is non-trivial in advance. However, it is 
possible to establish a statistical formalism of the space 
partition in a gas without an explicit description of the 
volume V. The application of the theory to the ideal gas 
allows us to identify the available volume in this simple 
system. This solution will provide a reference for under- 
standing more complex systems composed by interacting 
particles. 



III. GAS DESCRIPTION 

We consider a gas composed of N particles in a D- 
dimensional container of size V. The region y is a seg- 
ment, a surface or a volume for _D = 1, 2 or 3, respec- 
tively, although it will be called generically a volume. At 
low densities the state of the gas can be described using 
one-particle states. Since we are specially interested in 
the configurational state of the gas, the translational de- 
scription is here omitted and its well-known contribution 
to the fluid free-energy will be added later. 

The physical space occupied by the gas at a given in- 
stant is distributed among the N particles. To each parti- 
cle a certain volume v is assigned (so-called the available 
volume, see of the total space V. The statistical 
study of space partition can be performed based only on 
the constrains of normalization (there is a volume v per 
particle) and space-filling (the sum of all of them equals 
the total gas volume), without an exact specification of 
V in advance. We represent the configurational state of 
the gas by an occupation number distribution 7V„, where 
Nydv is the number of particles with available volumes 
between v and v + dv. The physically allowed distribu- 
tions must satisfy the following conditions 

[ N^dv, (1) 
Jo 

V = [ vNydv. (2) 
Jo 

Configurational entropy. The configurational entropy 
of the gas is represented by the Boltzmann-Planck ex- 
pression Sconf = kin W, where k is Boltzmann's con- 
stant and W the number of configurational microstates 
of the gas which reproduce a specific state Ny. Follow- 
ing the usual procedure, we subdivide the total volume 
into small cells which contain fractions dw of neighbor- 
ing states. In order to evaluate dw we turn the attention 
to well-known results of statistical mechanics. There are 
dl — Vd^p/h^ possible states for a particle confined in 
the volume V having a momentum p within a momen- 
tum volume d^p (h is Planck's constant). The number 
of states that corresponds to a physical volume d^x is 

d^x 

state number in {d^x, d^p) ~ dl . (3) 
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This result suggests that the fraction of configurational 
states of a particle having an available volume between 
V and V + dv (and any momentum) is proportional to dv 
with a proportionality factor V^^. 

Thus, a cell has a fraction of states dv/V and contains 
Ny dv particles for a given occupation set Ny . It is readily 
shown by well-known methods in statistical mechanics 
that 



(dv/V) 



N„ dv 



Nl 



n 



(dv) 



dv 



(4) 



The asterisk means that we must multiply over groups 
of single particle states instead of all states, according to 
the cell division. In the last term on r.h.s. of Eq. Q), V 
was removed from the product taking into account that 
J2l Nydv = N. We stress that the assignation of a vol- 
ume V to each particle does not imply a knowledge of the 
particle positions in the physical space or that the parti- 
cles should be at specific sites. In this sense, the present 
configurational description of a system of particles can 
be applied to different matter structures (gases, liquids 
and solids). 

The calculation of Sconf is straightforward. Using Stir- 
ling's theorem in the form ml « (m/e)™, we obtain 



S. 



conf 



= -k 



Ny\n 



NyV 

N 



dv . 



(5) 



where we have returned to an enumeration over all one- 
particle states instead of groups of states and replaced a 
sum over v by an integral. 

The Helmholtz energy. The free energy of the gas is 
the sum of three terms 



F = Ft, 



F, 



conf 



(6) 



Fconf — —TSconf IS thc configuratioual free-energy (be- 
ing T the temperature of the gas) and Ftrans the trans- 
lational free-energy, which in equilibrium is given by, 



^ trans — d I 



(7) 



where (3 = l/{kT) and A = (27rft2/3/r7i)i/2 ig the thermal 
de Broglie wavelength {h = /i/(27r) and m is the mass of 
one particle) . Uint is the sum of all energy contributions 
originated from the interactions among the particles. It 
is calculated as follows. 

Interaction energy. The average of the sum of two- 
particle interactions may be expressed as 



N r 

Uint = 2Y J ' 



where 



[V] 



Uy{uj) gvi(^) du) . 



(8) 



(9) 







^\ (0 


/[V] 



FIG. 2: Volumes involved in the evaluation of interaction 
energy between a reference particle (central point) with avail- 
able volume V (schematically represented by a circle) and the 
rest of the gas contained in V (shaded area). The volume cu 
centred in the chosen particle is used to evaluated the inter- 
action of this with the rest of the gas. [V] is the minimum 
volume ij which encloses the gas. 



Here Uv{oj) is the pair interaction potential and gv(w) 
the pair distribution function (p.d.f.) corresponding to 
particles with available volume v. Both functions are 
expressed as a function of the volume uj enclosed by the 
interaction distance r. 

Eqs. inj arise from considering the interactions 
between a particle having an available volume v and 
ngy{uj)dLu particles {n = N/V) located between the sur- 
faces of spherical volumes uj and to -\- duj centred on the 
chosen particle. The integration over oj takes into ac- 
count all possible interactions of this particle with the 
remaining ones in V . This evaluation is repeated (in- 
tegration over v) for considering the interaction of each 
particle with the remaining gas. A factor of 1/2 avoids 
the double count of binary interactions. The upper limit 
\V] of integration in Eq. © depends on the position of 
the reference particle within the gas volume and must be 
chosen as such that it completely covers V (see Fig. (SJ. 
However, it is assumed in Eq. that 4>y is independent 
from the particle position in V [since Uy{ijj) — > rapidly 
as w — > cxD, for typical fluids]. We must adopt the con- 
vention that portions of lj outside the volume V do not 
contribute to the integral. This prescription is applied to 
all integrals over ld of the present paper. Since ngy{uj) 
integrated over V gives iV — 1 (i.e., the total number of 
particles around a chosen particle with available volume 
w), the p.d.f. (7i,(w) must verify 



gy{ijj)dbj — V 



(10) 



The dependence of the pair particle potential with the 
available volume assumed in Eq. is not superfluous. 
It is well-known that many-body effects can modify the 
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interaction between two particles. For example, iV-body 
effects commonly soften the repulsive part of the two- 
body interactions in atomic and molecular fluids. In liq- 
uid hydrogen, the repulsion between two H2 molecule 
is modified when a third molecule (or more) is placed 
nearby, because this induces changes in the electronic 
clouds of the particles [s^l • A softening effect on pair po- 
tentials is also observed in liquid helium at high density 
|33l |. The confinement effects in the forces experimented 
by the particles can be implicitly included in effective 
pair potentials (or so-called pseudopotentials) [s^, for 
example with the adoption of a density-dependent pair 
potential ^3 ■ Alternatively, many-body effects in binary 
interactions could be represented in the present frame- 
work with a dependence of the interaction potential (uy) 
of a particle on the free space around it. It is worth 
noting that if the pair potential is independent from v, 
the interaction energy {4>y) of a particle can be even a 
function of the available volume throught the p.d.f. (gv)- 
When the pair potential is independent from the avail- 
able volume of the reference particle, Eqs. H8I9|I are re- 
duced to 

Ar2 AV] 

Uint = — I u{uj) g{uj) duj , (11) 



2V 



where 



The equilibrium of the particle distribution is deter- 
mined by the usual method. For low density gases, the 
minimization of F with the subsidiary conditions [Eqs. 
((TOl ] yields 



1. iV.A^ 



2V 



N 



a + 



dv 



-0, 
(14) 

where a and 7 are the Lagrange's parameters associated 
with conservations of the particle number and the vol- 
ume. Since Eq. p4|l must be verified for arbitrary func- 
tional variations {(5iV„}, we conclude that the equilibrium 
density is 



N4>v , Urnt 



2V 



N 



(15) 



Substitution of this expression in Eq. H13|l gives the equi- 
librium Helmholtz energy, which can be easily written in 
the form of a Euler-like equation 



F 



-7^ + aN- Ui, 



(16) 



The parameters 7 and a are closely related to the pres- 
sure P and the chemical potential /x. In order to demon- 
strate it, we first derive two useful identities. By intro- 
ducing Eq. ((T3)l in Eq. (P) we find that. 



g{uj) = {g,{w))^- g,{uj)N,dv. (12) 

g{ijj)duj is the probability to find a particle at a distance 
between r and r -\- dr (ld being the size of a ZJ-sphere 
of radius r) from another fixed particle regardless of its 
available volume. Therefore, this p.d.f., written in terms 
of the interparticle distance r, coincides with the radial 
distribution function g(r) typ ically used in theories of 
liquids and plasmas (e.g., |33)- Consequently, Eq. itTT)) 
is the standard expression of the interaction energy in 
a fluid. Therefore, the conventional form of the total 
potential energy of a one-component, classical fluid is 
automaticaly recovered from Eqs. JHEJ, as well as the 
p.d.f. g{r) is recovered from Eq. H12fl using the more 
detailed p.d.f. 5t,(w). 



IV. EQUILIBRIUM STATES 

The three contributions to the r.h.s. of Eq. JHJl may 
be joined together into a unique expression. The total 
free-energy of the gas in kinetic equilibrium now reads 



F = 1 








Jo 




V N J 


2V _ 



dv . 



(13) 



At low densities, high-order correlations in the particle 
distribution may be neglected so that ^1,(0;) and (there- 
fore) <j)y are independent from N^. The exact form of 
f/^(w) depends on the considered gas. We shall return to 
this point later (Sees. IVl and IVl|l . 



Adv 



(17) 



where A is the exponential term of Ny in Eq. (|15|l . Be- 
sides, since the current independent variables are the 
temperature, the particle number and the volume, we 
formally have a = a{T,N,V) and 7 = -f{T,N,V). Dif- 
ferentiating expression H17|l with respect to N (using the 
Leibniz's theorem [3^) we obtain 



dMCt 



N 



V 



Adv 



vAdv^ — 



(f>yAdv = . 



(18) 



Integrals in Eq. H18I) may be written in terms of A^, y 
and Uint ■ A little algebra leads to the first identity. 



NdNa - VdNl - dNUint = 



(19) 



The second one follows from the differentiation of Eq. 
(|17|l with respect to V . The result is 



Ndva ~ Vdvi - dvU.nt + ^ 
N 

"W 9v<f>vN,,dv = . 



(20) 



Here a term proportional to Ny^y has been dropped 
since it tends to zero in the thermodynamic limit (T- 
limit, N,V ^ 00 keeping N/V constant). Typically, the 
interaction potential Uy does not depend on V, so the 
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integral term in Eq. H2U|) may also be neglected in the 
T-limit. Then, the second identity takes the form of 



Ndva - Vdvl - dvU.nt + 



V 







(21) 



From the usual thermodynamic relations and Eq. 116(1 
it follows that 



^ = OnF — a + NdNce — Vd^j — d^Ui, 



(22) 



P = ~dvF = -f + Vdv7-Ndva + dvU,nt- (23) 
With the help of identities ((T!!|l and ((2U, we find that, 
H = a, (24) 



P = 7 + 



V 



(25) 



Therefore, a is identified with the chemical potential, 
while 7 equals to the pressure minus the density of inter- 
action energy. Substituting Eqs. H24|) and (|25|) into Eq. 
(|TB|l yields 



F = -PV + fiN 



(26) 



according to the Euler's relation p37j . 

The program for the application of the present the- 
ory is as follows. Given the pair potential of a fluid, the 
evaluation of Eqs. and 10) with the equilibrium dis- 
tribution [Eq. (|15|l ]. yields the chemical potential ^ and 
7 in terms of T, N, and V. Substitution of these so- 
lutions and Uint [which is a function of N and V, see 
Eq. ((SJ] into Eq. (|16(l gives the characteristic function 
F = F{T,N,V). As it is well-known, the availability 
of the thermodynamic potential F in terms of the inde- 
pendent variables T, N and V , provides us with the full 
thermodynamics description of the system. 



V. PERFECT GAS 

For a gas composed of non-interacting particles 
Uv{uj) = (jiv = Uint = 0. The chemical potential /i, 7 
and the space distribution are determined by substitut- 
ing Eq. H15|l into the conditions given in Eqs. and 
Hence we obtain, 



fcTln 



A^/j7 
1 - e-/3-r^ 



1 - (1 + l3-iV)e-'^'^^ 



Pie 



1 



(27) 



(28) 



(29) 



Since there is no interaction energy in the gas, P = 7. 
In the T-limit, equations of state (|27|l and 1(281) recover 
their usual forms {(3^ — > n) 



= kT In (nX^) , 
P = nkT 



(30) 



while the spatial distribution of particles is described by 



V 



exp {—nv) 



(31) 



Henceforth, Uy is called double density because its units. 
From Eqs. ((26() and l(3U() . the free energy takes the ex- 
pected form F = NkT [in (?iA^) - l] . 



A. Identification of v 

We have noted in i^that the elucidation of the avail- 
able volume per particle is non-trivial and difficult to 
find previous to the application of the equilibrium the- 
ory. The present application to a perfect gas shows that 
the volume V is distributed among the particles following 
a decreasing exponential law [Eqs. ((29() and ((31() ]. Most 
particles have available volumes going to zero, while few 
particles have a larger self space. This behavior is the 
same for all dimensions D — 1, 2 and 3. 

It is instructive to compare this result with the size dis- 
tribution of Voronoi regions associated to systems of non- 
interacting particles. The size distribution of Voronoi 
tesselations in one dimension is given exactly by a F- 
distribution [2^ 

2 

n^/"-"^"^) = ^(cnw)(^-i) eM-cnv) , (32) 
r(c) 

where c = 2D and F denotes the usual gamma function. 
For higher dimensions {D > 1) there are no exact results, 
although the size distributions can be even fitted by Eq. 
((32() for two and three dimensions ■ As we mentioned 
in Sec. El the Voronoi regions consist of all points closest 
to the particle chosen and, therefore, represent segments, 
poligons and polyhedra for D = 1, 2 and 3, respectively. 
Althought the Voronoi space partition could be consid- 
ered a good cantidate for the equilibrium distribution 
derived in Sec. |lH[Eq. by comparison of Eqs. 1(3 1() 

and l(32() . it is clear that the available volume per parti- 
cle adopted by the thermostatistical theory applied to 
perfect gases is not a Voronoi region. 

However, there is a simple interpretation of v for a 
perfect gas, as we now observe. We notice that the dis- 
tribution Uy obtained for a system of non-interacting par- 
ticles is equivalent (after a normalization) to the Hertz or 
nearest- neighbor distribution function H(r) associated to 
such systems [s^. In fact, the number riydv of particles 
with available volume between v and v + dv given by Eq. 
((31() is just the number H(r)dr of particles which have the 
center of the nearest particle lying at a distance between 



FIG. 3: Typical configuration of uncorrelated points in two 
dimensions. Sfiaded circles denote the available volumes of 
the objects (see text). 

r and r + dr, i.e., between the surfaces of D-dimensional 
spheres of sizes v{r) = v and v{r + dr) ~ v + dv centred 
on the reference particle [see Eq. (|38|l ]. This analogy 
suggests that the volume w of a particle in a perfect gas 
in equilibrium can be directly related to its distance to 
the nearest neighbor. Specifically, for a system of non- 
interacting bodies we propose the following identification 
of the available volume per particle, 

V = _D— dimensional sphere of radius given by the 
distance between the centers of the reference 
particle and its nearest neighbor. 

(33) 

Fig. O illustrates the available volumes associated to a 
two-dimentional collection of Poisson points. One notices 
that the proposal (|33|l does not establish an ownership 
connection between points of the space and each parti- 
cle, so that the space is not divided into distinct, non- 
overlapping regions. Nevertheless, we have not found any 
conflict between this identification of v and the theory. 
On the contrary, this proposal for v contains two useful 
advantages. First, as we shall show below, the present 
identification of v leads to a simple and accurate repre- 
sentation of the p.d.f. gy. Second, v gives an adequated 
notion of the free-particle space surrounding each particle 
and provides an useful tool for evaluating the perturba- 
tions experimented by atoms and molecules in diluted 
real gases. We will return to this point in Will 

Fig. 2] shows a direct test of Eq. (|31|l and the identifi- 
cation H33|l by computer simulation. Systems of iV = 10^ 
uncorrelated points were simulated in unitary spaces of 
one, two and three dimensions using a standard algo- 
rithm [33 . The exponential decreasing behavior of the 
size distribution of available volumes is clearly repro- 
duced by the numerical experiments. As a reference, one 



12 3 

Cell size 

FIG. 4: Results obtained from a simulation with 5 x 10^ Pois- 
sonian points in one, two and three dimensions. Exponential 
smooth curves correspond to the theoretical prediction [Eq. 
I3II1 I. The Voronoi' statistics given by Eq. 13211 is shown by 
comparison (dashed lines). 

can observe in the figure the notably different laws fol- 
lowed by the Voronoi's space partitions (dashed lines). 

B. Pair distributions 

From the identification H33() . the evaluation of the p.d.f. 
91,(0;) associated to particles with available volume v is 
immediate. We notice that v gives explicitly the distance 
to the first neighbor, which has a contribution gi^^ to gy. 
Since gi^\Lo)dLO represents the probability of finding the 
nearest neighbor between the surfaces of spheres w and 
Lu + duj, one find that n gi^\Lj')duj' must be zero for 
uj < V and one otherwise. Therefore, 

gip(,o) = -S{u;-v), (34) 
n 

where 5{uj — v) is the Dirac delta function. On the other 
hand, v does not give any information about the positions 
of the remaining particles in the gas except that they are 
outside of a spherical volume v centred on the reference 
particle. Besides, these neighbors are uncorrelated with 
the reference particle. Thus, their total contribution to 
gv must be one for lu > v and zero in the contrary case, 
which constitutes a Heaviside step function Q{uj — v). 
Therefore, the p.d.f. of reference particles with available 
volume V can be expressed in the form 

gv{uj) — —S{uj — v) + Q{ll) — v) . (35) 
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Fig. [51 confirms our identification of It is also 

worthwhile to notice that Eq. H35|l reproduces the ex- 
pected value of the averaged p.d.f. [Eq. ((T^ ] for a set of 
non-interacting particles, 



(7„(a;)ne-™dz; = e-"'^ + l-e-"'^ = 1, (36) 



which expresses that the positions of the particles are 
completly uncorrelated and, therefore, the probability 
density of finding neighbors (without any specification 
of the available volume of the reference particle) is uni- 
form and equals unity. We notice also that g^^'>{uj) = 
{gi^\uj)) = e~"'^ represents the first neighbor distribu- 
tion for random configurations of non-interacting points. 
The nearest-neighbor distribution function H{r), which 
was introduced in W Al can be calculated from the rela- 
tion 



H{r) = n— — g'- > [uj{r)] . 
ar 



(37) 



Its evaluation reproduces the well-known results of a per- 
fect gas j4^]. 



Hir) 



duj{r) 
n — - — exp I 
dr 



;(r)] 



(38) 



For = 3, it is a classical result derived by Hertz |4ll |. 

In closing this section, it is worth stressing that all ter- 
modynamics and spatial statistical, well-known results of 
classical perfect gases in one, two and three dimensions, 
are self-consistently derived from the theory developed 
in Sees. IlIII and IIVI The identification of the available 
volume given in (|33|l plays a central role in the present 
application of the theory, and so we believe that it should 
be regarded as the most plausible interpretation of v for 
systems of non-interacting particles. 




1 2 3 4 5 

Normalized radial distance 

FIG. 5: Pair distribution function for particles with avail- 
able volume V — (1.1 ± 0.1)7i~^ in a two-dimension ideal gas 
as a function of the dimensionless distance x = r^/n (units 
of the mean distance between particles). The peack around 
X = 0.59 is the contribution of the nearest neighbor. The 
results are for A'^ = 10000 particles randomly distributed in 
a unitary square surface (662 particles were selected with the 
chosen v). The insert graph shows An = Joidv — l)irxdx, 
i.e., the difference between the actual particle number within 
a radius x with respect to that one corresponding to a (micro- 
scopical) uniform particle density. The unitary jump (within 
the numerical errors) in a; = 0.59 confirms that the peak of 
Qv corresponds to one particle (the nearest neighbor). 



VI. GASES OF HARD ELASTIC PARTICLES 

Systems formed by hard bodies constitute simple fiuid 
models in which impenetrable particles interact solely via 
hard-core repulsions. These models have played a major 
role in the liquid state theory, as the behavior of dense 
fluids is dominated by the exclude-volume effects asso- 
ciated with the repulsive forces of its constituents. The 
main advantage of these models is the simplicity of the 
pair potential, which can be defined by 

where a is the size of a Z?-sphere of radius d, d being the 
diameter of a particle. The aim of this section is to apply 
the theory to hard-particle systems, which can provide a 
guidance for the treatment of more realistic fluids. We 
concentrate on the limit of low densities, although the 
results derived for D = 1 will be valid at all densities. 
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FIG. 6: Schematic representation of the available volume 
u of a hard rod (central shaded rectangle) in presence of its 
nearest neighbor. The volume Vi is the ("ideal") available 
volume of the central body if the bodies are fully penetrable 
(non-interacting), according to the interpretation 1331 . 



As a consequence of the hard-core interaction, the iden- 
tification H33I) of the available volume of a particle is not 
valid for the present fiuid. In fact, we notice that the 
space occupied by a particle is inaccesible to the other 
ones, therefore this region must be a part of the avail- 
able volume of the particle not shared with the others. 
Consequently, the space occupied by neighbor particles 
yields a reduction of the available volume of a particle 
with respect to that (vi) calculated for fully penetrable 
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TABLE I; Quantities used in the study of D- dimensional 
hard-body fluids. io{r), size of a _D-sphere of radius r; a, 
volume of a hard-body of diameter d; rj, packing fraction; a, 
exclusion sphere of the repulsive interaction; a* , reduction of 
the available volume per hard-body; b = a — a* (see text). 
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FIG. 7: Same as Fig. El but for D — 2. The exclusion region 
a (circle closed by dotted points) is a forbidden region to the 
particle center of neighbors. The available volume of a hard 
particle (w) is lower than the perfect-gas value (vi) (see text). 



(non-interacting) objects, so that v < Vi. As an ansatz, 
we propose that the available volume of a hard _D-sphere 
is given by 

w = - a* , (40) 

where a* is the reduction of the available volume of a 
hard particle with respect to the perfect-gas value 
As we shall show later, the value of a* can be evaluated 
using the virial theorem. 

On the other hand, one expects that the p.d.f. of 
a perfect gas [Eq. H35|l ] is a good approximation to the 
p.d.f. of non-ideal fluids at very low densities. If we take 
into account the relation H4U|I . the p.d.f. of a diluted 
hard-body fluid can be written as 

g^{oj) ^ -S{uj - Vi) + e{uj ~ v^) . (41) 
n 

From Eqs. ©, H39|l and (|41|1 . we find the mean potential 
energy of a particle with available volume v, 



oo, 
0, 



V <b 

V > b 



where 



(42) 



(43) 



From Eq. (|15|) then results that A^^, is zero ioT v < b and 
has an exponential dependence on the available volume 
for V > b. Clearly Uint = and then P = 7. At the 
T-limit, evaluations of Eqs. (|1I2|I yield 



kT 



nb 



1 - nb 



In 



1-nb 



nkT 
l-n6' 



= Q{v — b) 



1-nb, 





n{v — b) 


j cxp 


1-nb 



(44) 



(45) 



(46) 



The free energy is easily derived from Eqs. H16(l . I|44() and 
and reads 

F = NkT{lTi[inX")/{l-nb)]-l}. (47) 

As we anticipated, the virial theorem may provide the 
value of b and, thus, the elucidation of the factor a* at low 
densities. According to Eq. (|45|l , b represents the second 
virial coefRcient in the equation of state P — nkT{l + 
bn -\- ...) (since ^ — P) and, therefore, its value is well- 
known for hard D-spheres. If we denote with a the size 
of a hard-particle, then b = a, 2a and 4cr, for _D = 1, 
2 and 3 |42|. Therefore, from Eq. (|43|l and the known 
values of a, it follows that a* — b — a/2. Explicit values 
of interest are given in Table 

The averaged nearest-neighbor p.d.f. can be calculated 
substituting the first term on the r.h.s. in Eq. H41|l and 
n^ given by Eq. H4t)|) into Eq. H12|) . Thus, we obtain 
that, 



,(^)(.)^e(.-a)(^) 



exp 



n{Lj — a) 



1 — nb 



(48) 



Further, the p.d.f. for all neighbors except the nearest 
one can be derived in similar way from the second term 
on the r.h.s. of Eq. (|41|) . The result is 



<?('')(w)==e(cj-a)h-exp 



n{ijj — a) 
1 — nb 



(49) 



The total averaged p.d.f. is composed by the sum of these 
two contributions. 

Knowledge of ^(^^(cj) permits us to evaluate the 
nearest-neighbor distribution function H(r) from Eq. 
(|37|l . For convenience, we introduce the dimensionless 
distance x = r/d and the packing fraction 77 = na, which 
gives the fraction of the total volume occupied by the 
hard bodies. Then H{x) = d H{r), and we obtain 



H{x) 



H{x) 



H{x) 



1-77 



cxp 



-2-q [x - 1) 
1-^ 



877 a; 
l-2r; 



24?7x2 
1 -4r; 



exp 



exp 




x > 1 (50) 



X > 1 (51) 



x>l (52) 



for _D = 1, 2 and 3, respectively {H = at a; < 1). 
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FIG. 8: Dimensionless double density of a hard-rod fluid for 
values of the packing fraction rj — na — 0.1, 0.5, and 0.8. 
Dotted curves correspond to the theoretical prediction [Eq. 
I46l l1. and solid lines represent the results obtained from nu- 
merical simulations with 2 x 10* hard particles. The Voronoi' 
statistics is shown by comparison (dashed lines). 



A. Hard-rod fluid 

The properties of a fluid of hard rods {D = 1) are 
known exactly, hence this model offers an appropriate 
test for our evaluations. Thus, it can be verified that the 
free energy given in Eq. (|47|l is rigorously exact for D = 1 
and, therefore, all the equations of state are correct. In 
particular, Eq. gives the exact equation of state for 
hard rods as derived earlier by Tonks [43l |. 

As we have seen in the case of a perfect gas, the present 
theory not only gives the thermodynamic properties of 
the fluid, but also provides information on its microstruc- 
ture throught the doble density n„. Fig. |S1 shows riy for 
several values of the packing fraction = 0.1, 0.5 and 
0.8, as obtained from Eq. Due to the effect of the 

core repulsions [see Eq. no particle may have an 

available volume lower than 6, so that the distribution 
Uy is characterized by an abrupt fall in v = b {nv = 77 in 
the figure). The distribution n„ becomes more and more 
peaked as the gas concentration increases. This predic- 
tion is confirmed by numerical simulations (Fig. (SJ. 

It may surprise us that the p.d.f. given by Eq. H41() 
is only an approximated expression (its form was adopted 
from the perfect-gas solution) , while its application leads 
to correct results as previously described. This is due 
to the fact that the approximation made on gy rests on 
the second term on the r.h.s. in Eq. H41|l (it ignores pair 
correlations between a particle and its second, third, etc.. 



neighbors), which is not required for the determination 
of the double density riy. In a linear gas only the p.d.f. of 
the nearest neighbor is relevant, being correctly given by 
the first term on the r.h.s. in Eq. l|IT|l . Thus, Eq. (gSJ 
is exact for D = 1 and leads to the correct distribution 
H{r) of the nearest neighbor, as given in Eq. H50|l . This 
result was derived earlier by MacDonald |4J| . 

Finally, we recall that the present theory yields exact 
results for the thermodynamic and geometric structure of 
hard rods. It gives an important support to the identifl- 
cations 133|l and 140|l of the available volume per particle. 



B. Hard disks and spheres 

Not a single exact result is known for D > 1, but a 
number of numerical solutions and analytical approxi- 
mations are available [4^. Inspection of the free energy 
[Eq. (|17|l ] and equations of state (0^ and (|^ evaluated 
with values of b given in Table reveals that our results 
arc found to give the correct thermodynamical properties 
of hard discs and spheres at low densities. 

In order to verify the double density given in Eq. 
(|46|l . we study the nearest neighbor distribution func- 
tion H{x) which is derived from it. For this purpose, 
Eqs. (|51|l and (|52|) are contrasted with the analytical 
expressions derived by Torquato , which show an ex- 
celent agreement with Monte Carlo simulation data i^. 
We use Torquato's expressions valid for densities lower 
than the freezing density (the freezing packing fraction 
is ry/ « 0.69 and 77/ « 0.49 for discs and spheres, re- 
spectively). Figs. 1^ and [TUl show that Eqs. H51() and 
(|52|l give the correct asymptotic behavior of H{x) at low 
densities, and capture the salient features of H(x) for 
moderate densities. The agreement between these results 
and Torquato's expressions is within 1% everywhere at 
packing fraction rj < 0.01 and 77 < 0.002, for 13 = 2 and 
D = 3, respectively. 

In closing, the low-density, asymptotic results of hard 
discs and spheres are correctly reproduced by the present 
theory. We find again a support for the identification of 
the available volume per particle, closely related with the 
distance to the nearest neighbor. 



C. Van der Waals fluid 

In the preceding subsections, we have studied the ap- 
plication of the theory to systems composed by particles 
with infinite hard cores. We shall consider now the pre- 
vious model but with the inclusion of an attractive long- 
range contribution in the potential. Because attractive 
forces do not restrict the accessible space to a particle, 
the identification (|40|l of the available volume per par- 
ticle (with a* given in Table D) should remain unaltered 
for this model. Instead of trying to figure out the pair 
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potential itself, one may directly define as 



FIG. 9: Nearest neighbor function H{x) of a liard-disc fluid at 
packing fraction 77 = 0.02, 0.1, and 0.2, in terms of the reduced 
distance x = r/d (d, diameter of a particle). Solid curves 
result from Eq. 151L dashed lines correspond to accurate fits 
up to freezing {j] < 0.69) 40], and dotted lines are results for 
fully penetrable particles. 
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FIG. 10: Nearest neighbor function H{x) of a hard-sphere 
fluid at packing fraction r) = 0.01, 0.02, and 0.05, in terms 
of the reduce distance x = r/d {d, diameter of a particle). 
Solid curves result from Eq. I52II . dashed lines correspond to 
accurate fits up to freezing (ry < 0.49) 40], and dotted lines 
are results for fully penetrable particles. 



00, 
-e. 



V <b 

V > b 



(53) 



where e is a constant parameter (with units of energy 
times volume) and b is given by Eq. H43|) . Eq. H53|) is 
essential for obtaining the Van der Waals equation and is 
the simplest choice to test the relation (I25|l between the 
pressure and the parameter 7 when Umt is non-zero. 

In fact, this form for is based on the suitable short- 
range repulsion and an attractive interaction whose range 
is finite but long compared to the interparticle spacing. 
An appropriate attractive potential with these character- 
istics is the general Kac potential 



,Ko 



(54) 



with Jp°° J{uj)duj — 1. Here 7* represents the Dth power 
of the usual inverse-range parameter. The evaluation of 
the attractive-interaction energy per particle is immedi- 
ate. From Eqs. iflT)) and one obtains (in the 
T-limit) that, 



,Ka 



= l*Jh*{v + a*) 

n 



e j J{lu)lli . (55) 



In the limit 7* 0, 0^'°^ — e and we recover the at- 
tractive term in Eq. I|53|l . This result shows that the 
attractive energy of a particle results from its simultane- 
ous interaction with a large number of particles, where 
the contribution of the nearest neighbors is negligible re- 
gardless their distances to the reference one. As a conse- 
quence, (j>y°''^ is equivalent for all particles independently 
from their available volumes. 

The interaction energy given in Eq. (|53ll . after sub- 
stitution in Eq. (|15|l and normalization, reproduces the 
double density given by Eq. I|46|) . In the T-limit, Eqs. 
(P) and 10) lead to 



fi — kT 



nb 



l-nb 



In 



nX 



D 



1-nb 



en , 



(56) 



while expression (|45|l is obtained for 7. The evaluation of 
Eq. (jSJ yields a non-zero interaction energy with density 
Uint/y = — en^/2. Besides, from previous results and 
Eq. (|16|l . the Helmholtz energy takes the form of 

F ^ NkT {In [71X^/(1 - nb)] - l} - eNn/2 . (57) 

From differentation of F with respect to V one obtains 
the well-known Van der Waals equation of state, 



P 



nkT 
1 — nb 



en 



(58) 



Relation (|25|) is thus confirmed beyond the trivial case 
(P = 7) obtained for perfect and hard-body fluids. The 
Van der Waals model constitutes a special case where at- 
tractive forces have not influenced the value of 7, which, 
however, is strongly governed by the inflnitely repulsive 
interactions. 
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D. Order in the fluid structure 

We notice that the double density offers a way to mea- 
sure the order of the packing of hard particles. In fact, 
we can define the fluctuation of the available volume as 
the standard deviation 

where (...) represents an average in the ensemble [earlier 
applied in Eq. (fT^ J , namely, 

1 

ifv) = - fvnydv , (60) 
Jo 

fy being a function of the available volume per particle. 
Obviously, (v) — n^^ for all fluids [see Eq. ©J. Accord- 
ing to Eq. (|46|l for hard-body fluids, one obtains 

Aw ~ {v) — b , (hard bodies) . (61) 

For comparison, the space fluctuation of classical perfect 
gases is as large as the mean volume per particle for any 
temperature and density [c.f. Eq. (|^ ]. 

Av — (v) , (perfect gas) . (62) 

This points out the irregular geometric structure that 
characterizes random collections of non-interacting ob- 
jects. 

The volume fluctuation Av given by Eq. H61|l . summa- 
rizes the behavior of riy with the density for hard rods. 
Large deviations of the available volumes from the mean 
value occur at low densities, and the hard-rod fluid has 
the inhomogeneity of a perfect gas in the limit — > 
{Av ^ n^^). As expected, the fluctuations decrease at 
high densities and vanish in the limit of close-packing 
rj = 1. Eq. (|61l) is valid for discs and spheres at low 
densities. 

Classical methods for determining the order in an 
isotropic packing of hard D-spheres are based on geomet- 
rical and topological techniques. The former rests on the 
identification of peaks in the radial function g{r), which is 
found to be unsatisfactory because, among other reasons, 
it is difficult to determine when the peak appears |5lj |. 
Topological methods use the so-called order or topolog- 
ical parameters, such as like-coordination numbers |49j 
and the number of faces in a cell division 

[llE3- They 

commonly require an appreciable amount of calculation 
divided in two steps, first, an evaluation of configura- 
tions of the system in equilibrium states (usually based 
on numerical simulations), and, then, the execution of 
an algorithm which evaluates the cell division (typically 
a Voronoi tesselation) . 

In constrast with the previous methods, the present 
formalism offers, eventually, a direct and precise way to 
characterize the order in fluid structures in terms of the 
fluctuation of the available volume per particle. Thus, for 
example, Av = (v) represents the high density fluctua- 
tions present in very diluted gases. On the contrary, uni- 
formily ordered configurations should be characterized by 



space fluctuations noticeably lower than (v). Although 
the application of Av is here limited to hard particles 
at low densities, additional investigation could extend its 
use to real fluids. 



VII. DILUTED HYDROGEN PLASMA 

Finally, as the main example of this work we consider 
a low density gas formed by hydrogen atoms in an elec- 
trically neutral background of free electrons and protons. 
For simplicity, we consider a strongly ionized equilibrium 
state where the atom-atom interaction can be ignored, 
and only the charge effects upon the atoms are taken 
into account. Each atom is considered as a polarizable 
particle and its available volume is approximated by the 
perfect-gas value Our aim focuses on the descrip- 

tion of the neutral component of the gas. We emphasize 
that the results obtained from this model are qualitative 
and illustrative only. A detailed study of a mixture of 
gases in equilibrium is in progress and its results will be 
presented elsewhere. 

The equilibrium density given in Eq. H15|l can be easily 
generalized to particles with internal states following the 
procedure of Sees. IIIII and IIVI The double density of 
atoms in the energy level j and with available volume v 
[hereafter, atoms (j, v)] is given by 

Tl 

nj,v = ^ 9j exp [-(3 (e^ + Ae-jy + jv - /i)] , (63) 

where gj and ej are the multiplicity and eigenenergy of 
the level j for an isolated atom, n is the total density, and 
A and fi the thermal wavelength and chemical potential of 
the atomic component, respectively. The factor Aejy can 
be interpreted as the energy change of a particle (j, v) due 
to interactions with the charges. Specifically, we obtain 

Aejy = i {ne(t>jv,e + np(t>jy^p) , (64) 

with 

4>jv,c^ / Uji,,c(t^) 5ju,c(w) , (65) 
Jo 

where the subscript c refers to either electrons or protons, 
Ue and Up are the number densities of free electrons and 
protons, Ujy^c is the pair potential of the interaction of an 
atom (j, f) with a charge c, and gjy^di^) represents the 
probability density of finding a charge c in the surface 
of a spherical volume oj centred on an atom {j,v). We 
have omitted in Eq. (|64|) the energy change of atoms in 
the level j due to their perturbations over the charges, 
because it is only a small correction (proportional to the 
abundance of atoms) in Aejy for the physical conditions 
assumed. 

At large separations, charges induce a dipole moment 
in the hydrogen atom. The pair potential associated to 
polarization of an atom by either an electron or a pro- 
ton can be represented by the screened Debye potential 
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obtained by Redmer and Ropke [5^ (in the form used in 



Upoi{r) 



1 + r/rs \ It It, 



I] + r2 



(66) 



where e is the elemental charge, the screening length, 
Cij the polarizability of an atom at the level j, and Ij = 
asj[(7j2 + 5)/4]^/2 (with j the main quantum number 
and as the Bohr radius) its mean radius. 

At short separation between an atom and a proton the 
Coulomb repulsion between protons prevails |53| | 



uc{r) = — 1 + - 
r V a/ 



-2T/a 



(67) 



This potential includes a plasma screening function with 
a characteristic length a — ■^UBlj/h- 

Due to Pauli's exclusion principle, there is a strong 
repulsion when the electronic wavefunctions of two atoms 
begin to overlap and this is often expressed as 



ui2{r) 



Cl2 
„12 



(68) 



We use this potential to simulate the interaction be- 
tween a free electron and an atom. The coefficient C12 
is harder to evaluate. From Hindmarsh et al. [s^l we 
adopt C12 = 9 X 10~^^^j erg cm^^, where we consider 
only the wavefunction size (of characteristic radius Ij) of 
the bounded electron. 

On the average, roughly half of atoms have an electron 
as its nearest neighbor and the other half a proton. The 
fraction of atoms having other atoms as nearest neighbors 
is negligible in a highly ionized material. Therefore, the 
total density of atoms can be written as 



E 

3 



njjdv 



(69) 



where n^^"* and n'j^J are the double densities of atoms 
which have a proton or an electron as their nearest neigh- 
bor. Both densities are given by Eq. (times 1/2), 
although they differ in the energy term Aejy. For the 
first and second atom sets we have, 



Ae 



jv 



[neUe{uj) + npMp(cj)] du) 



[rieUeiuj) + rip Up (w)] duj 



(70) 



(71) 

respectively, where Up — uc + Upoi and Ue = M12 + "po/ 
are the total atom-proton and atom-electron potentials, 
respectively, and r = (Sv/Att)^/^. In these evaluations, 
Eq. iT?^ was adopted for the species (electrons or pro- 
tons) which let the nearest neighbor of an atom {j,v), 
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FIG. 11: Double density (above) and energy (below) curves of 
the first ten levels of hydrogen atoms (solid lines from left to 
right), which have an electron as the nearest neirghbor. The 
results are shown as a function of the available volume per 
atom, and correspond to T = 15000 K and mass density p = 
10"'' gr cm~^. The dotted line shows the energy continuum 



edge for electrons and protons. 



while we approximate gjv^ci^^) = 0(lu — v) as the p.d.f. 
associated to the other perturber class. 

From Eq. (|53|) and the division of atoms into two 
groups (according to the type of nearest perturber), the 
internal partition function of the atomic hydrogen is 
given by 



-(3'yv 



-/3Ae ' 



4- e 



dv . 



(72) 

This partition function includes a state sum over the 
physical space as a direct consequence of the configu- 
rational parameter v used in the gas description. When 



the interactions are negligible (Ae^,„' 



Ae/,^ = 0), the 



present partition function reduces to Z = Z^/n where 
Zq = gje~^'^^ is the conventional form of the internal 
partition function of a perfect gas. 

Populations and effective internal energies of atoms as 
functions of the available volume per particle are calcu- 
lated for a gas with a temperature T — 20000 K and 
mass density p = 10~^ gr cm~'^. This is representative of 
conditions found in the atmospheres of lukewarm white 
dwarfs (effective temperature around Te// = 15000 K). 
Non-ideal effects on equations of state are negligible and, 
particularly, (3^ = n = 10^^ °^ cm~'^. A standard evalua- 
tion of the chemical equilibrium yields Ug, ^ Up = lO^^ ''^ 
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cm ^ and nn = 10^^-^^ cm ^. 

The effective internal energy of an atom is composed 
by the sum of its electronic energy ej (that chosen of 
an unperturbed atom) and the interaction energy Aej^ 
which condenses the plasma effects. Lower panels in Figs, 
im and El show the effective internal energy as a func- 
tion of its available volume for the lowest atomic levels. 
Atoms with great volumes are, in practice, unperturbed 
and their internal energies are very close to those of an 
isolated atom. As v decreases the atoms become progres- 
sively more bound due to the attractive long-range inter- 
action with charges. This interaction has a stronger v- 
dependence for high-lying levels due to a greater shearing 
of the electronic wavefunction, and accordingly, a greater 
polarizability. 

The Coulomb repulsions and Pauli's exclusion effects 
preclude atoms from having very low available volumes. 
These short-range interactions increase the effective in- 
ternal energy of atoms and bound states successively 
merge into the continuum with decreasing the available 
volume. Differences between energy curves of Figs. El 
and El at low v are due to the fact that Coulomb re- 
pulsion between protons [Eq. (|67|l ] is appreciably softer 
than the electron-atom interaction [Eq. Ht)8|) ]. In par- 
ticular, atoms with an electron as the nearest neighbor 
have an abrupt increasing of interaction energy at low 
available volumes. Notice that linear Stark effects are 
not taken into account because we use energy eigenvalue 
(cj) of unperturbed particles. However, since the linear 
Stark effects only spread over the levels without shifting 
their centers of gravity, the main effect on the internal 
energy of the atoms is due to shifts caused by plasma 
polarization and short-range interactions considered in 

Density curves in Figs. II H and [T^ (upper panels) show 
that atoms exist only for available volumes within a cer- 
tain domain, which depends on the internal state j. The 
double densities take maximum values around the mini- 
mum of the internal energy curves. The absence of par- 
ticles with very low volumes v is due to the fact that 
atoms with high internal energies are not stable. Indeed, 
according to the energy curves, an atom with very low 
available volumes is so strongly perturbed that its elec- 
tron becomes free in the plasma. Atoms in high-lying 
bound states are more strongly affected by the presence 
of plasma electrons and ions, because their electrons are 
loosely bound. Since the Coulomb repulsion uc is softer 
than the potential ui2, atoms with a proton as the near- 
est neighbor may survive at lower available volume than 
those with an electron as the nearest neighbor. However, 
this difference does not have an influence on the total 
populations of each group of atoms. 

The population of atoms in all internal states declines 
at high V volume {v « 10~^^'^ cm~^) because actu- 
ally there are no available volumes per particle much 
greater than the mean volume per particle in the gas, 
(v) = n^^ = 10^^^ °^. Mathematically, this is a con- 
sequence of the term exp(— /?7ti) in the double density 




Log v [cm~^] 

FIG. 12; Same as Fig. 1111 but for atoms with a proton as 
the nearest neighbor. 

expression [Eq. H63|l ]. Such factor plays the role of a 
probability of occupation of available volumes. Large 
volumes v represent single-particle, configurational states 
with smaller occupational probability than the little ones. 

Fig. ll3lshows total densities (i.e., w-averaged) of atoms 
in differents internal states, together with predictions ob- 
tained for non-interacting particles. Low-energy bound 
states have essentianly unperturbed populations. On the 
contrary, atoms in high energy levels are so strongly per- 
turbed that their equilibrium populations decrease with 
respect to a hypothetical gas of unperturbed atoms. Un- 
der current conditions, states with main quantum num- 
bers larger than « 30 are not populated because the 
atom, in order to survive, requires a volume v much larger 
than the mean volume per particle and the probability of 
existence of such space becomes very small. Notice the 
smooth dependence of the non-ideal effects on atomic 
density with the quantum main number. This contrasts 
with evaluations of internal partition function based on 
an (often used) abrupt cutoff at certain upper level. The 
present evaluation of internal partition function contains 
useful advantages derived from the free-energy minimiza- 
tion technique: (i) convergence of the internal partition 
function, (ii) incorporation of particle-interaction models 
in a consistent statistical mechanical way. 

A variety of gas models have been developed to com- 
pute the equilibrium state of astrophysical fluids, and 
which include the influence of non-ideal gas effects upon 
the internal partition function of bound states [ll ls^lsH . 
However, to our knowledge, this is the first time that 
combined and self-consistent evaluations of atomic popu- 
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FIG. 13: Populations of atomic levels as a function of the 
quantum main number for T — 20000 K and logp = —7. 
Open circles show results for atoms with an electron as the 
nearest neighbor, while points correspond to atoms with a 
proton as the n.n. (actually both atom groups have very 
similar populations). The dashed line shows results based on 
an ideal gas (i.e., for Ae^^' = Aej^' = 0). 



lations and internal effects on bound states are performed 
differentiating groups of atoms under different plasma- 
perturbations. We consider this as an important advan- 
tage of the present formalism, which can be useful for 
more elaborated models. Notice also that the data de- 
rived from the population equation could be used 
to set up thermodynamic functions and to obtain optical 
quantities. This is not however within the scope of this 
paper. 



VIII. CONCLUSIONS 

Classical statistics of gases in equilibrium has been ex- 
tended taking into account a new one-particle variable 
which assigns a space v to each particle. We call this 
variable the available volume per particle since it roughly 
represents the free space around it. We develop the the- 
ory in terms of a Helmholtz free-energy model for one- 
component gases at low densities where the particle in- 
teractions are assumed to be binary. The equilibrium 
state is calculated from the energy minimization with re- 
spect to the occupation number Ny of u-states, subjected 
to particle number and volume constrains. The volume- 
closure introduces a Lagrange's multiplier, here called 
7, which regulates the space partition (i.e., iV„) and is 
equivalent to the pressure minus the density of interac- 
tion energy in the gas. A universal expression [Eq. (|15|l ] 



gives the occupation number Ny for dilute gases at equi- 
librium. When the pair potential is specified, Ny contains 
both thermodynamics and structure information of the 
fluid. We have illustrated the application of the theory 
to the perfect gas and to fluids with inflnitely repulsive 
interactions, for one, two and three dimensions. All ther- 
modynamic and geometric-structure properties such as 
the nearest neihgbor distribution function of these sys- 
tems have been correctly reproduced. For hard bodies 
in two and three dimensions, the results are valid at low 
densities. It was also shown that the Van der Waals 
model can be derived from general Kac potentials in a 
straighforward way. 

The identification of v has been perfomed from the 
application of the theory to specific cases. For nonin- 
teracting particles the available volume is equivalent to 
the spherical volume enclosed by the distance between 
the centers of the particle and its nearest neighbor. In 
hard-body fiuids, the core-repulsions reduce the value of 
V with respect to that for fully penetrable particles. At 
the limit of low densities, the explicit form of v for hard 
bodies is derived from the virial theorem. Hard particles 
provide a guidance for the treatment of soft-cores present 
in realistic fluids. In these systems, we expect that the 
allowed values of v will be a compromise between values 
associated with noninteracting particles and those corre- 
sponding to hard particles. Additional work should be 
done on this topic. On the other hand, there is no ev- 
idence that the identification of v can be modified by 
attractive interactions. This is reasonable because at- 
tractive forces do not restrict the accessible space to a 
particle. 

For non-interacting particles Ny follows an exponen- 
tially decreasing law, such that most particles (a fraction 
1 — e^^ « 0.63) have spaces v lower than the mean vol- 
ume per particle (l/f^)• In conexion with the structure of 
fluids, we have showed that the double density Uy could 
be a useful tool to characterize the order in fluids, par- 
ticularly, through the standard deviation of the available 
volume V. 

The main application of the formalism concerns with 
the atomic component of a partially ionized hydrogen 
gas. We found a statistically coherent and detailed evalu- 
ation of populations and internal energy changes of atoms 
as a function of the available volume per particle. Pre- 
liminary calculations show that the distinction of atoms 
under different plasma-perturbations emerges naturally 
from the theory. This finding provides sufficient confi- 
dence for further extension of the formalism to consider 
a complete description of gaseous mixtures, where non- 
ideal effects on optical properties may be investigated. 
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